Contact resistance dependence of crossed Andreev reflection 
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We show experimentally that in nanometer scaled superconductor/normal metal hybrid devices 
and in a small window ol contact resistances, crossed Andreev reflection (CAR) can dominate 
the nonlocal transport for all energies below the superconducting gap. Besides crossed Andreev 
reflection, elastic cotunneling (EC) and nonlocal charge imbalance can be identified as competing 
subgap transport mechanisms in temperature dependent four-terminal nonlocal measurements. We 
demonstrate a systematic change of the nonlocal resistance vs. bias characteristics with increasing 
contact resistances, which can be varied in the fabrication process. For samples with higher contact 
resistances, CAR is weakened relative to EC in the midgap regime, possibly due to dynamical 
Coulomb blockade. Gaining control of crossed Andreev reflection is an important step towards the 
realization of a solid state entangler. 



Quantum mechanically entangled pairs of particles are 
a major building block of quantum computation and in- 
formation processing. A natural source of entangled elec- 
trons is a BCS-type superconductor where the Cooper 
pairs (CP) form spin singlet states. The two electrons 
of a Cooper pair can be spatially separated into two dif- 
ferent metallic leads in a nonlocal process called crossed 
Andreev reflection (CAR) 0, [3, i, 3 . 

At temperatures (T) well below the superconducting 
transition temperature, T c , and for bias potentials be- 
low the superconducting energy gap A, electrons from a 
normal metal contact (N) can enter the superconductor 
(S) only as Cooper pairs by a process known as Andreev 
reflection (AR). In this local process a hole is reflected 
into the same N to conserve momentum. If two normal 
metal contacts, Nl and N2, are spatially separated by 
less than the coherence length £, the two electrons form- 
ing a CP can originate from different normal contacts, 
see Fig. QJa). This process opens an additional nonlo- 
cal conduction path known as CAR. An inverse process 
was proposed as the basis of a solid-state entangler: the 
electrons of a CP are split between the two leads while 
retaining their entanglement from the superconductor. 
However, this method of creating entangled particles can 
be accompanied by two additional processes that lead to 
correlated signals on Nl and N2, but not to entangle- 
ment. In the first, a single electron from Nl can reach 
the other contact N2 by elastic cotunneling (EC) 
see Fig. [TJb). In the second, called nonlocal charge im- 
balance (CI), electrical charge can be transferred to the 
second contact by the diffusion of quasi-particles gener- 
ated by finite temperatures or finite bias. 

Recently, the relative strength of these subgap pro- 
cesses was the subject of extensive theoretical work. 
Standard BCS theory predicts that to lowest order in the 
tunneling rates, CAR and EC exactly cancel in normal 
metal/insulator/superconductor (NIS) systems at low T 
and bias This cancelation can be lifted for higher 



transmissions Q, by spin-active interfaces Q or ferro- 
magnetic contacts J7 ] , by disorder, or by electron-electron 
interactions I 111 . Il2j . It has been suggested that 
Coulomb interaction and the electromagnetic environ- 
ment are crucial for CAR to be the dominant process 

Experimentally, first signatures of nonlocal transport 
processes have been reported only recently on multi- 
terminal hybrid structures. We can group the experi- 
ments into two categories: while EC and CI dominate 
the subgap transport in devices with highly transpar- 
ent contacts [3, [HI , CAR is more prominent in samples 
with lower transparencies, depending also on the applied 
bias 



mm 



These individual findings suggest that the 
contact resistance plays a significant role in determining 
which process is dominant. 

Here we present four-terminal nonlocal experiments on 
a series of samples of similar geometry but with NIS junc- 
tions of different resistances and with non-magnetic con- 
tacts to avoid stray field effects. We demonstrate that 
in a small window of contact resistances CAR can domi- 
nate the nonlocal transport for all subgap bias voltages at 
low temperatures, and that CAR is suppressed for larger 
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FIG. 1: (Color online) Schematics of (a) crossed Andreev 
reflection and (b) elastic cotunneling. (c) SEM image of a 
typical sample. A current / is applied between normal contact 
Nl and the superconducting contact SI, while the nonlocal 
voltage U n \ is detected between N2 and S2. 
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barrier resistances. 

In our experiments we perform four-terminal nonlocal 
measurements: a bias U between the normal metal con- 
tact Nl (injector) and the grounded superconductor lead 
SI causes a current / as shown in Fig.QJc). We measure 
both, / and the nonlocal potential difference U n \ between 
the normal contact N2 (detector) and the second super- 
conducting contact S2. In a CAR process a particle is 
injected from Nl and one of opposite charge is generated 
in N2. In Fig. QJa) such a process is shown schemati- 
cally for U < 0, where electrons are injected from Nl 
and holes are created in N2. For each U a potential of 
opposite sign builds up on N2 to ensure the zero current 
boundary condition of a voltage measurement. In EC 
and CI processes, however, the charges on the injector 
and detector have the same sign and U n \ has the same 
polarity as the bias. In practice, we use ac modulation 
techniques and investigate the nonlocal differential resis- 
tance R n \ = dll n i/dl « U^/I^, where 'ac' denotes the 
amplitude of the corresponding variable. i? n i is negative 
for CAR and positive for all other subgap processes. 

The samples are planar multi-terminal hybrid struc- 
tures, see Fig. [He). They were fabricated on thermally 
oxidized Si wafers by e-beam lithography and using the 
angle-evaporation technique in ultra-high vacuum with a 
base pressure of 10 _10 mbar. Because of the low spin- 
orbit interaction, we chose Al as superconducting mate- 
rial. The Al layer has a thickness of 50 nm and a width 
of 200-300 nm. We oxidized the Al between 3 and 15 min 
in situ in an oxygen atmosphere with a pressure in the 
range of 0.1 to 12 mbar to obtain the required tunnel bar- 
riers. Without breaking the vacuum, we tilted the sample 
and deposited Pd as normal metal contacts with a typ- 
ical thickness of 30-40 nm and a width of 90 nm. In our 
Al/AlOz/Pd-samples we obtained edge-to-edge distances 
d between the N electrodes of ~ 50 — 220 nm. 

For each device the characteristic energy A* denotes 
the bias potential of maximum conductance of the injec- 
tor contact in superconducting tunneling spectroscopy 
(STS) measurements, see Table I. We note that A* de- 
pends on the barrier characteristics and does not coin- 
cide with the bulk energy gap of A « 200 //eV. We mea- 
sured a critical temperature T c « 1.2 K and a resistivity 
p A i = 1.2 - 2.6 ^Ocm of the Al layers at T = 4.2 K. 
We find an elastic mean free path t e \ sa 20 — 50 nm and 
£ w 180 — 270 nm at cryogenic temperatures [2l[. There- 
fore the Al is in the 'dirty limit' of superconductivity. 

We performed nonlocal ac measurements using lock-in 
amplifiers at a frequency of ~ 10 Hz in a He 3 -cryostat 
with a base temperature of 0.23 K. The dc bias applied 
to Nl, [/dc, is modulated by C/f c ps 12^Vrms. Since the 
junction resistances vary strongly with the applied bias, 
leakage currents and capacitive cross-talk can produce 
spurious signals which might resemble the expected char- 
acteristics. To exclude cross-talk we repeated all mea- 
surements at various frequencies. For all samples shown 



here we find that the in-phase part of U n \ is indepen- 
dent of the measurement frequency, while the capacitive 
part Y n \ is negligible, see Fig. [21(c). For comparison we 
simultaneously measure U n \ in a DC configuration with 
10s delay between points (not shown). The dc coincide 
with the ac results, but suffer from a much worse signal- 
to-noise ratio. Another source of spurious ?7 n i signals is 
resistive leakage into the detector. We reduce this leakage 
by using voltage amplifiers with 1 GQ input impedance 
and check for each device that a change to 100 MO has 
no effect on our results. We note that for T > T c all 
nonlocal signals vanish, implying that they are related to 
the superconductor and not to the measurement setup or 
inhomogeneous current paths. 

In the following we present experiments on three sam- 
ples A, B, and C, which have a similar geometry, but 
differ in the contact resistances. The characteristics of 
each sample are summarized in Table [U 

Figure [21(a) shows R n \ as a function of the bias voltage 
Ud c at various temperatures for sample A. At base tem- 
perature we observe a negative nonlocal resistance i?„i for 
all subgap voltages with R n \ ss — 1.0 fi. This value does 
not change with temperature up to T rs 0.5 K. When 
increasing T further R n \ increases and becomes positive 
for T > 1.05 K. For voltages larger than A*/e, R n \ is 
positive at all temperatures and increases with T. At 
1.6 K, i.e. for T > T c , i? nl » for all biases. The con- 
tact characterization of sample A at base temperature 
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FIG. 2: (Color online) Results of sample A. (a) R n i vs. bias for 
various temperatures with Nl as injector and N2 as detector, 
(b) STS characterization of contacts Nl and N2 at 0.23 K. 
The dashed lines label the bias corresponding to A* of the 
injector, (c) Frequency dependence of the in-phase nonlocal 
signal Uni and the capacitive part Y n \. 
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TABLE I: Properties of samples A. B and C. t e \: elastic mean 
free path of the Al strip at T = 4.2 K; £: superconduct- 
ing coherence length; Ri, R d : injector and detector barrier 
resistance at T = 1.5 K; A*: characteristic energy of the 
injector; d: injector-detector distance; -RiAi, RdAd- resis- 
tance area product of the injector and detector; asymmetry 
a = (R d A d - RiAi)/{R d A d + RiAi). 
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is shown in Fig. [DJb) and exhibits a significant subgap 
conductance for both, Nl and N2. The dashed lines rep- 
resent the bias corresponding to the energy A* of the 
injector, which coincides with the sign change of R n \. 
Within the energy gap both barriers exhibit nearly the 
same normalized conductance values and both contacts 
have similar resistances above T c , 0.4 kfi and 0.5 kfi. Fig- 
ure [2fc) shows the frequency dependence of the nonlocal 
signal as discussed above. 

In Fig. [3] the bias dependence of i? nl of sample B is 
plotted for several temperatures. The shape and ampli- 
tudes of the curves are very similar to those of sample A 
in Fig. HJa), as is the temperature dependence. However, 
in sample B at base temperature we observe a local max- 
imum for Udc ~ with R n \ w 0, whereas for finite subgap 
biases we find R n i < with two minima of R n i ~ — 1.1 fi 
at Udc ~ ±0.235 mV. Above T c the nonlocal signals van- 
ish for all bias voltages. The inset shows the STS charac- 
terization of the individual contacts at 0.23 K. In contrast 
to sample A, the normalized injector and detector con- 
ductances differ considerably in the gap and are smaller 
than in sample A = 0.8 kfi and Rd = 5.5 kfi). 
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FIG. 3: (Color online) Results of sample B. R n \ vs. bias 
at various T with Nl as injector and N2 as detector. The 
inset shows the STS characterization of Nl and N2 at 0.23 K. 
The dashed lines indicate the bias corresponding to A* of the 
injector. 



Bias dependent measurements for sample C are shown 
in Fig. [UJa). At base temperature a prominent local 
maximum in the middle of the superconducting gap de- 
velops with a positive maximum of R n \ « 30 fi. For 
0.13 mV < \U dc \ < A*/e with A* w 0.31 meV, i?„i is 
negative with minima at ±0.23 mV with R n \ w — 13 fi. 
The sign changes at Udc ~ ±0.13 mV, independently of 
temperature. For |C7dc| > A*/e the signal is positive and 
tends to zero for increasing bias, in contrast to samples A 
and B. With increasing temperature all nonlocal signals 
tend to zero, independently of the bias, with R n \ w for 
T > T c . At zero bias the signal decreases monotonically 
with increasing T and vanishes already for T > 0.5 K, 
see Fig. HJc). A finite nonlocal signal only develops for 
T well below T c . From the STS measurements shown 
in Fig. 0Jb) we find that the two contact resistances are 
quite different and much larger than in the previous sam- 
ples (R = 2.6 kfi and R d = 11.2 kfi). 

Based on the discussion above we interpret R n \ < 
in the superconducting energy gap as CAR being the 
dominant process, while R n \ > is attributed to EC or 
CI. The latter two processes can be distinguished by the 
temperature dependence of R n \ at a given bias. EC and 
CAR are reduced with increasing temperature [l]| Gil , 
while CI increases strongly up to T c due to the decreased 
superconducting energy gap and the broadened thermal 
distribution, and vanishes for T > T c [HI]. At zero bias 
and low enough temperatures CI can be neglected. 

For sample A we conclude that CAR is the dominant 
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FIG. 4: (Color online) Results of sample C. (a) 7? n i vs. bias 
for various T with Nl as injector and N2 as detector. The 
solid lines are guides to the eye and the dashed lines show 
the voltage corresponding to A* of the injector, (b) STS- 
measurements of contact Nl and N2 at 0.23 K. (c) R n i at zero 
bias as a function of T. 
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nonlocal transport process at base temperature and at all 
subgap biases. From the characteristic strong increase of 
R n \ to positive values with T, we conclude that nonlo- 
cal CI becomes important for T > 0.5 K. These findings 
also apply for sample B, except around zero bias, where 
CAR and EC approximately compensate each other at 
low temperatures. At finite subgap bias, EC is weakened 
as reported in [l6| and the characteristics become simi- 
lar to sample A. In sample C, EC not only compensates 
CAR, but even dominates the nonlocal midgap transport. 
As in A and B, CAR dominates at finite subgap bias. In 
all samples CI dominates for bias potentials above the 
energy gap. In contrast to the first two samples, CI is 
strongly suppressed in sample C for \Udc\ >> A*/e. 

Our experimental results on different samples allow us 
to examine which device parameter determines the R n \ 
characteristics. We now discuss the following possible 
parameters: 1) the distance d between injector and de- 
tector 2) the resistance of the injector and detector con- 
tacts 3) the resistance-area product, RA, of the contacts, 
and 4) the RA asymmetry a between injector and de- 
tector. First we note that d does not change the char- 
acteristic shape of the R n \ curve [l6j |. In addition, it is 
well established, both theoretically and experimentally, 
that all nonlocal subgap signals at low bias are reduced 
monotonically with d 0, 14, [HI, [Til]. Since the junc- 
tion separation is largest in sample C, which exhibits 
the largest nonlocal signals, we conclude that as long 
as d < £, the characteristics of the nonlocal response is 
determined by another parameter. Our results suggest 
that larger resistances of the injector and detector con- 

\ compa 

0,113 

curves consistent with our data when we use RA as the 
relevant parameter. This is remarkable since in [l6| the 
geometry and materials are very different. In contrast to 
RA, we find no correlation between the RA asymmetry 
a = (RdAd — RiAi) I (RdAd + R\Ai) and the systematic 
evolution of our data from sample A to sample C. From 
this discussion we conclude that the shape of the R n \ 
curves of our samples is determined by RA. 

Our qualitative explanation for this dependence is 
based on a model in which the two normal contacts are 
electromagnetically coupled, which can favor CAR over 
EC and lead to R n i < 0, depending on the symme- 
try of the coupling fill ]. This model requires dynami- 
cal Coulomb blockade (DCB) on the two NS junctions 
in union with the complex impedance of the S contacts. 
This leads to a blockade on the superconductor which 
suppresses CAR, because the net charge — 2e of a CP 
is added to the superconductor and the corresponding 
charging energy has to be supplied. Since CAR also van- 
ishes for completely transparent contacts this leaves 
only a small window of transparencies where CAR can be 
observed at zero bias. In contrast to CAR, no net charge 
is added for EC, which is therefore essentially unaffected 



tacts suppress CAR and favor EC. A comparison with 
similar experiments in the literature 

EE Ea yields i? n i 



by the blockade and dominates the subgap transport for 
low transparencies around zero bias. 

DCB depends strongly on the environmental 
impedance and thus on the resistance of the tun- 
neling contacts, Rt [HI- It was shown that DCB can 
occur also for Rt < h/2e 2 « 12kf2, and that decreasing 
Rt leads to a reduction of the DCB effects RA 
is the natural parameter to characterize a finite area 
tunnel junction with many parallel channels since it 
scales with the average channel resistance. A large 
number of weakly transmitting channels may lead to a 
small Rt while retaining DCB for each channel (lit ]. 

The systematic change of the measured nonlocal resis- 
tance vs. bias characteristics in our samples can thus be 
understood as follows: the contacts of sample A have the 
lowest resistances and the blockade is lifted, which allows 
CI and CAR to develop. The latter is favored over EC 
possibly due to an electromagnetic coupling of the injec- 
tor and detector characteristic for all our samples. For 
sample B with slightly larger contact resistance, DCB 
sets in and weakens the dominance of CAR around zero 
bias which leads to a reduction of the negative nonlocal 
resistance. At a finite subgap bias EC is weaker because 
the process can not be elastic anymore. With resistances 
of a few kf2 as in sample C, DCB becomes strong and 
inhibits CAR relative to EC. We note that (local) AR 
is suppressed stronger in the DCB regime than CAR on 
two independent (but spatially close) contacts [2(| , which 
might explain the larger R n \ values in sample C. Since for 
CI a net charge is transferred to the superconductor, the 
DCB picture is also consistent with the suppression of CI 
in sample C at large bias. 

In conclusion we find in our experiments that for small 
RA CAR can dominate the transport characteristics for 
all subgap biases, whereas above the superconducting en- 
ergy gap the CI rates are important. We have shown that 
CAR is suppressed for larger RA around zero bias. We 
qualitatively explain the various subgap transport char- 
acteristics by different dependences of the relevant pro- 
cesses on dynamical Coulomb blockade, which depends 
crucially on the RA. 
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